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Andrew J. Baker 1 , Richard I. Davies 1 , Matthew D. Lehnert 1 , Niranjan A. Thatte 1 ' 2 , William D. Vacca 1 , 
Olivier R. Hainaut 3 , Matt J. Jarvis 2 ' 4 , George K. Miley 4 , and Huub J. A. Rottgering 4 

1 Max-Planck-Institut fiir extraterrestrische Physik 
Postfach 1312, D-85741 Garching, Germany 

e-mail: {ajb, davies .mlehnert , thatte , vacca}@mpe .mpg . de 

2 University of Oxford Astrophysics 

Keble Road, Oxford OX1 3RH, United Kingdom 
e-mail: m j j Oastr o . ox . ac . uk 

3 European Southern Observatory 

Alonso de Cordova 3107, Casilla 19001, Vitacura, Santiago, Chile 
e-mail: ohainautOeso . org 

4 Sterrewacht Leiden 

Postbus 9513, 2300 RA Leiden, the Netherlands 
e-mail: {miley , rottgeri}@strw . leidenuniv . nl 

Received 7 March 2003 / Accepted 2 June 2003 

Abstract. As a prerequisite for cosmological studies using adaptive optics techniques, we have begun to identify 
and characterize faint sources in the vicinity of bright stars at high Galactic latitudes. The initial phase of this work 
has been a program of K s imaging conducted with SOFI at the ESO NTT. From observations of 42 southern fields 
evenly divided between the spring and autumn skies, we have identified 391 additional stars and 1589 galaxies 
lying at separations A8 < 60" from candidate guide stars in the magnitude range 9.0 < R < 12.4. When analyzed 
as a "discrete deep field" with 131 arcmin 2 area, our dataset gives galaxy number counts that agree with those 
derived previously over the range 16 < K a < 20.5. This consistency indicates that in the aggregate, our fields 
should be suitable for future statistical studies. We provide our source catalogue as a resource for users of large 
telescopes in the southern hemisphere. 
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1. Introduction ing role and subtle interactions of key physical processes 

such as feedback, mass loss, merging, and secular evo- 
Understanding the mechanisms that drive galaxy forma- among otherg Disentangling these pr0C esses in 
tion and evolution is a central objective of modern as- nearby systems ig &lready extremely difficu i t . Attempts to 
trophysics. While the structures and dynamics of nearby extend ^ effort tQ newly _ formed galaxies require data 
galaxies constitute a valuable fossil record, these are much whoge Hgh gpatial resolution and dynamic range can only 
less informative than the properties of galaxies observed be achieved by 8m _ dass telescopes using adaptive optics 
at the epochs of their formation and early evolution. By (AQ) techniques to work at or near their diffract ion limits, 
studying galaxies' structures, dynamics, fundamental scal- 
ing laws, and distributions and rates of star formation AO relies on the detection of, and compensation for. 
over a range in redshift, it is possible to trace the evolv- the distortions of wave f r onts from a bright pointlike ob- 

~~Z ; ~ , . _ , iect. For science targets that are themselves too faint to 

bend offprint requests to: A. J. Baker, ajb@mpe.mpg.de ' . . 

* n , , . , . permit waveiront sensing, or too dmuse to sustain cm- 

Based on observations obtained at the European ° 

Southern Observatory, Chile, for programmes 66.A-0361 cient wavefront correction (see Rousset D2H1, a bright 

and 68.A-0440. Note: the entirety of Table El is only ( V < 13 )> nearby (AO < 30") natural guide star (NGS) is 

available in electronic form at the CDS via anony- required. Exactly how bright and nearby depends on the 

mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via local atmospheric conditions; the above values correspond 

http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A/ to the typical regime in which the wavefront sensor is not 
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photon-limited and the off-axis Strehl ratio at 2.2 /im is 
not too severely reduced (e.g., Le Louarn et al. I1998[l . 
In the absence of a nearby NGS, higher-order wavefront 
corrections using a laser guide star (LGS) can suffice, pro- 
vided a moderately bright (V ~ 16) reference star lies 
relatively close (AO < 90") to the science target for deter- 
mination of the lowest-order tip-tilt correction. Switching 
from LGS to NGS mode trades a gain in sky coverage 
for a loss in the highest possible Strehl ratio, however, 
since the LGS samples atmospheric turbulence in a coni- 
cal (rather than cylindrical) volume. The phase error due 
to such focus anisoplanatism increases with mirror diam- 
eter as £> 5 / 6 (e.g., Le Louarn et al. 119981 and references 
therein). Until the advent of multiple-LGS systems that 
can correct for this effect (e.g., Tallon & Fov ll99Uf) . AO 
systems on 8m-class telescopes will therefore reach the 
highest Strehl ratios only when operating in NGS mode. 

In order to conduct cosmological studies with the 
present generation of AO technology (e.g., Larkin et al. 
2000] Davies et al. IM)T1 Glassman et al. 1217)21 Davies et 
al. [2003 ) , it is necessary to identify distant galaxies in 
the vicinity of bright guide stars. Unfortunately, most ex- 
isting surveys either avoid bright stars (e.g., the Hubble 
Deep Field), or like DENIS and 2MASS (Epchtein et al. 
1997 Skrutskie et al. H997[) are too shallow to detect very 
many targets at cosmological distances. We have there- 
fore begun a program to characterize the extragalactic 
sources lying close to bright stars at high Galactic lati- 
tudes. Similar work underway by other authors (Larkin 
& Glassman 1999 ) has focused on fields easily observable 
from the Keck telescopes; our targets are specifically cho- 
sen to be at southern declinations suitable for observations 
with the NAOS-CONICA (Rousset et al. IT9"9"8l Lenzen et 
al. IT9"9"8"|l and SINFONI (Thatte et al. IT9"9"g|> instruments 
on the ESO Very Large Telescope (VLT). 

The initial phase of our program has been a campaign 
of (non-AO-assisted) near-IR imaging, since information 
about source magnitudes and surface brightness profiles 
will be most useful at the wavelength of eventual AO op- 
eration. For deep imaging studies in particular, the avail- 
ability of seeing-limited data for comparison will allow the 
empirical determination of the surface brightness selection 
effects influencing diffraction-limited data (e.g., whether 
a given faint "point" source is only the nucleus of an 
unremarkable extended galaxy). Blind diffraction- limited 
imaging of random fields, in contrast, remains vulnera- 
ble to unknown biases and uncertainties in the surface 
densities of true point sources. We have already followed 
up our near-IR imaging with optical imaging and multi- 
object spectroscopy of many of our fields; discussion of 
these data is deferred to future papers. 

2. Observations and data reduction 

2.1. Target selection 

Tabled lists the identifications, coordinates, and magni- 
tudes of the 42 bright stars that define our near-IR imag- 



ing sample; Fig^shows their distribution on the sky. The 
majority of these targets were selected from the USNO- 
A2.0 catalogue (Monet et al. I1998|l according to crite- 
ria designed to ensure optimal performance of a near-IR 
science instrument assisted by a visible- wavelength (e.g., 

0. 45 - 1.0 fim for NAOS-CONICA) wavefront sensor: 

1. magnitude 10.3 < R < 12.4, i.e., bright enough in the 
visible for wavefront sensing to be robust; 

2. color B — R < 1.1, i.e., bluer than that of a type Gl 
star, in order to maximize the number of photons on 
the wavefront sensor while minimizing the amount of 
near-IR light scattered onto the science detector; 

3. Galactic latitude \b\ > 15°, in order to minimize ex- 
tinction and contamination by large numbers of fore- 
ground stars; and 

4. declination in the range -44° < S < -13°, so that AO 
observations from Cerro Paranal (at latitude ~ —24°) 
can be obtained at low airmass. 

Because the USNO-A2.0 catalogue's star-galaxy separa- 
tion is not perfect (e.g., USNO-A2.0 0450-40377908 = 
ESO 436-G001), our search was not entirely automated. 
To check that favored objects were in fact stars, we used 
either the Tycho-2 catalog (H0g et al. I2000|l or the pub- 
lically available ESO Imaging Survey (EIS) data. Visual 
inspection of the J-band images for EIS Wide patches A 
(Nonino et al. ITM^ . B (Prandoni et al. ITM^ . and C and 
D (Benoist et al. 1999) allowed us to confirm that no ad- 
ditional bright stars or large foreground galaxies would 
block our view of more interesting background objects; 
we also tended to favor fields with one or more discern- 
able /-band sources lying within 30" of the central stars. 
In a small number of cases, we relaxed our magnitude 
and color criteria to include brighter and redder stars lo- 
cated within 20" of radio sources detected in the NRAO 
VLA Sky Survey (NVSS: Condon et al.\MB- Our goal in 
imaging these fields (noted in the last column of Table 0) 
was to take advantage of the fact that radio galaxies often 
trace overdensities in the large-scale structure, seemingly 
at all redshifts (e.g., Hill & LillyES]] Best 1277)771 Kurk et 
al.EHni Pentericci et al. 127)17)1 McLure & D union l2T7)T]l . 
Subsequent to our observations, the release of the USNO- 
B1.0 catalogue (Monet et al. 2003) allowed us to verify 
that only one of the stars in our sample (SBSF 05) has a 
proper motion as large as 50masyr~ 1 . 1 

2.2. Observations 

We obtained our observations using SOFI (Moorwood et 
al. HTTTsI on the ESO New Technology Telescope (NTT), 
over the course of six nights in February, March, and 
October 2001 (see Table 01. Conditions were photomet- 
ric for two of the three nights in each season, allowing us 
to obtain good photometry for all of our fields; K s seeing 
ranged from 0.5" to 1.0" with a median ~ 0.7". Our choice 

1 The preliminary photometry in the USNO-B1.0 catalogue 
shifts a few additional stars to B— R > 1.1, as TableQindicates. 
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Table 1. Bright star fields. Columns are (1) field identifier by which the sources in Tableware indexed (2) USNO- 
B1.0 designation for the bright star (Monet et al. lMI5|l (3) & (4) USNO-B1.0 J2000.0 coordinates of the bright star, 
accurate to ±0.2" (5) & (6) USNO-B1.0 proper motions of the bright star in masyr -1 , relative to the YS4.0 reference 
frame described by Monet et al. (2003) (7) & (8) USNO-B1.0 B and R magnitudes of the bright star, averaged over 
two epochs and accurate to ±0.3 mag (9) K s magnitude of the bright star, from our SOFI imaging (10) date(s) on 
which the field was observed with SOFI, as indexed in Tabled (11) total SOFI integration time, in minutes (12) note 
as to whether the field falls inside an EIS Wide patch or near an NVSS source. 



Field 


USNO-B1.0 


a (J2000.0) 


S (J2000.0) 


u (a) 


u(S) 


B 


R 


K s 


Night (s) 


At 


Note 






(2) 




(3) 








(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


SBSF 


01 


0703-0002281 


00 


H 


38.84 


-19 


37 


11.0 


+8 


-6 


12.1 


10.0 


7.80 


6 


30 


NVSS 


SBSF 


02 


0601-0007980 


00 


44 


31.88 


-29 


52 


30.3 


+24 


-8 


11.9 


11.1 


10.09 


6 


30 


EIS B 


SBSF 


03 


0600-0008265 


00 


-±'j 


20.62 


-29 


56 


46.0 


+32 


+6 


12.9 


12.1 


11.11 


6 


30 


EIS B 


SBSF 


04 


0604-0008038 


00 


45 


28.05 


-29 


31 


40.1 


-12 


+4 


13.2 


12.4 


11.33 


4 5 


20 


EIS B 


SBSF 


05 


0604-0008436 


00 


48 


11.23 


-29 


34 


14.2 


+92 


-60 


12.7 


11.6 


10.11 


6 


20 


EIS B 


SBSF 


06 


0605-0009459 


00 


~>0 


"}4 7fl 


-29 


26 


32.0 








12.6 


11.0 


9.90 


4 5 


30 


EIS B 


SBSF 


07 


0604-0008993 


00 


52 


02.66 


-29 


30 


48.1 


+24 


+40 


13.0 


11.7 


10.90 


4 5 


30 


EIS B 


SBSF 


08 


0605-0009658 


00 


52 


18.88 


-29 


27 


17.8 


-12 


-12 


13.1 


11.2 


10.81 


4 5 


30 


EIS B 


SBSF 


09 


0606-0009140 


00 


5*4 

oo 


1 7 9D 


-29 


22 


29.0 


-2 


-26 


12.2 


11.4 


10.21 


6 


30 


EIS B 


SBSF 


10 


0605-0009783 


00 


53 


17.95 


-29 


25 


22.6 


+26 


-42 


12.0 


11.1 


9.97 


5 


30 


EIS B 


SBSF 


11 


0656-0061729 


05 


36 


33.64 


-24 


19 


56.5 


+4 


-10 


11.8 


11.1 


10.08 


4 5 


30 


EIS C 


SBSF 


12 


0653-0066004 


05 




"}7 so 


-24 


36 


33.8 


+2 


-20 


12.1 


11.5 


10.26 


4 5 6 


30 


EIS C 


SBSF 


13 


0654-0065981 


05 


41 


49.93 


-24 


35 


32.2 








10.9 


10.8 


10.46 


5 6 


20 


EIS C 


SBSF 


14 


0767-0069418 


06 


07 


06.34 


-13 


13 


37.1 


+ 12 


-26 


9.8 


9.0 


7.99 


5 6 


20 


NVSS 


SBSF 


15 


0734-0203992 


08 


44 


00.22 


-16 


34 


01.1 


-4 





11.3 


11.8 


11.00 


1 


30 




SBSF 


16 


0705-0208711 


09 


24 


52.77 


-19 


26 


17.0 


-12 





11.0 


11.2 


10.89 


1 


30 




SBSF 


17 


0683-0253521 


09 


47 


44.79 


-21 


37 


12.7 


-8 


-12 


12.0 


11.6 


10.70 


1 


30 


EIS D 


SBSF 


18 


0682-0261875 


09 


40 


dp, QQ 


-21 


45 


13.3 








12.2 


11.4 


10.69 


1 


30 


EIS D 


SBSF 


19 


0682-0262596 


09 


51 


44.77 


-21 


43 


28.0 


-26 


+2 


13.2 


11.9 


11.06 


2 


30 


EIS D 


SBSF 


20 


0696-0235495 


09 


51 


50.64 


-20 


20 


14.8 


+8 


-34 


12.4 


12.4 


10.46 


2 


30 


EIS D 


SBSF 


21 


0686-0238319 


09 


-,o 


32.44 


-21 


21 


30.1 


-12 


+2 


12.7 


11.9 


11.33 


2 


30 


EIS D 


SBSF 


22 


0696-0236655 


09 


55 


14.78 


-20 


20 


03.3 


-12 





12.4 


10.6 


10.23 


1 


30 


EIS D 


SBSF 


23 


0528-0285326 


10 


09 


22.50 


-37 


11 


54.8 


-8 





11.1 


11.3 


11.27 


2 


30 




SBSF 


24 


0599-0250386 


10 


40 


26.20 


-30 


00 


36.5 


-14 


+8 


11.2 


11.6 


11.40 


2 


30 




SBSF 


25 


0577-0367919 


12 


24 


02.32 


-32 


14 


35.7 


-18 


-6 


10.3 


10.4 


10.28 


2 3 


30 




SBSF 


26 


0593-0288091 


12 


49 


55.31 


-30 


38 


18.7 


-26 


-2 


10.4 


10.4 


10.00 


2 3 


30 




SBSF 


27 


0582-0321519 


12 


55 


37.48 


-31 


42 


41.3 


-10 


-24 


11.4 


11.6 


10.97 


1 


30 




SBSF 


28 


0478-0346964 


12 


56 


14.14 


-42 


09 


10.9 


-8 


+4 


10.4 


10.5 


10.79 


2 3 


30 




SBSF 


29 


0580-0344020 


13 


03 


18.33 


-31 


54 


27.3 


-4 


-14 


10.8 


10.9 


10.90 


2 


30 




SBSF 


30 


0550 09S5RSQ 


13 


04 


37.76 


-34 


56 


30.8 


u 


1 U 


11.4 


1 1 7 


1 so 


I 






SBSF 


31 


0460-0297795 


13 


28 


32.45 


-43 


58 


03.4 


-14 


-2 


11.7 


12.1 


10.57 


1 


30 




SBSF 


32 


0562-0303721 


13 


28 


59.77 


-33 


42 


25.0 


-28 


+16 


10.3 


10.3 


10.05 


2 


30 




SBSF 


33 


0568-0379911 


13 


40 


31.54 


-33 


07 


52.6 


-12 


-16 


10.7 


10.8 


10.76 


2 


30 




SBSF 


34 


0582-0339867 


13 


46 


25.24 


-31 


45 


47.5 


+8 


-8 


11.4 


11.6 


10.97 


1 


30 




SBSF 


35 


0745-0822065 


21 


00 


25.58 


-15 


26 


12.4 


-14 


-12 


12.1 


10.4 


8.40 


5 6 


30 


NVSS 


SBSF 


36 


0615-0931505 


22 


14 


36.74 


-28 


25 


31.6 


-16 


-16 


11.7 


10.6 


9.36 


5 


30 


NVSS 


SBSF 


37 


0501-0829581 


22 


43 


04.43 


-39 


49 


29.3 


+28 


-22 


12.8 


11.8 


10.55 


6 


30 


EIS A 


SBSF 


38 


0498-0817047 


22 


47 


06.77 


-40 


10 


01.3 


+18 


-8 


11.5 


10.8 


9.29 


6 


30 


EIS A 


SBSF 


39 


0504-0833030 


22 


49 


34.23 


-39 


33 


05.3 


+18 


-36 


12.9 


11.1 


9.78 


5 


20 


EIS A 


SBSF 


40 


0501-0831035 


22 


49 


49.32 


-39 


53 


15.0 


-36 


-28 


12.4 


11.8 


10.64 


5 


30 


EIS A 


SBSF 


41 


0498-0817689 


22 


50 


21.28 


-40 


07 


38.6 


+18 





12.5 


11.2 


10.72 


5 


30 


EIS A 


SBSF 


42 


0714-0852904 


23 


29 


55.77 


-18 


35 


54.1 


+6 


+4 


11.0 


9.6 


7.90 


6 


30 


NVSS 



of observing strategy was guided by our desire to optimize 
the detection of faint sources at small angular separations 
from the bright star in each field. We used SOFI's small- 
field objective lens, giving 0.145" x 0.145" pixels and a 
2.47' x 2.47' field of view; larger pixels would only have 
hindered modelling of point spread function (PSF) wings 



close to the star, while a larger field of view would only 
have yielded additional sources too far from the star to be 
useful for future AO studies. In order to minimize mem- 
ory effects in the HgCdTe array, we used short detector 
integration times (2-4 s, depending on the magnitude of 
the star) and spent a total of one minute on source be- 
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Right Ascension (hours) 



Fig. 1. Positions of the 42 southern bright star fields observed in this survey, denoted by open circles. The solid line 
shows the Galactic plane; the dotted lines to either side indicate Galactic latitude |6| = 15°. The four EIS Wide patches 
from which some of the fields were selected are indicated by boxes (not drawn to scale). 



Table 2. Log of observations. 



Night 


Date 


Conditions 


(1) 


(2) 


(3) 


1 


2001 Feb 15 


photometric 


2 


2001 Feb 17 


photometric 


3 


2001 Mar 07 


thin cirrus 


4 


2001 Oct 06 


thin cirrus 


5 


2001 Oct 07 


photometric 


6 


2001 Oct 08 


photometric 



fore dithering the telescope to a new position. After imag- 
ing the field at ten such positions- each located randomly 
within a "jitter box" of dimensions 40" x 40"- we switched 
to a different bright star field for another ten-frame se- 
quence. We would return to a particular field again later 
in the same (or a different) night, in order to allow the 
diffraction spikes of the (Alt-Az mounted) NTT to rotate 
on the sky and thus expose any faint objects lying un- 
derneath. All of our fields were observed for a total of at 
least 20 minutes, and most were observed for a total of 30 
minutes (see Table 0). During each of the four photomet- 
ric nights, we interspersed multiple observations of six or 
seven faint standard stars (Persson et al. [1998) in order 
to determine the flux scale. 



2.3. Data reduction 

We reduced the data using custom scripts within the 
NOAO IRAF package fTodv ITMHj) . For each set of ten 
consecutive frames, we made an estimate of the (sky and 
dark current) background by masking out objects and- 
after further rejection of high and low pixels according to 
an estimated variance- taking a median. Following back- 
ground subtraction, each frame was divided by a dome 
flat field (the difference of lamp on and lamp off frames) 
that had been corrected for the dependence of the array 
shade pattern on incident flux. 2 Values for bad pixels were 
then interpolated. We next addressed the problem of low- 
level array cross talk, which resulted from the bright star 
in the field and manifested as excess counts in two verti- 
cal stripes. To make a more precise measurement of this 
rather weak effect, we coarsely aligned the ten frames, av- 
eraged them, and measured the median of each column 
in the mean image. The resulting values were then sub- 
tracted from the corresponding columns in the individual 
unaligned frames. Finally, we performed a fine alignment 
to compensate for irregular behavior by the image dero- 
tator during our runs. Although the derotator always as- 
sumed the correct position for the initial frame in each 
set of ten, it often remained at the same position during 
the subsequent nine. We therefore rotated and translated 
frames 2-10 in each set to align them with frame 1. For 

2 This type of "special" dome flat is described in Version 1.3 
of the SOFI User's Manual. 
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the final image combination (of all 20-30 frames taken for 
each field), after rejection of high and low pixels according 
to an estimated variance, we calculated the mean of the 
remaining pixels at each point. The angular resolution in 
the final images was calculated as the FWHM of a Moffat 
(1969) profile fit to the central star in each field. 

For each observing run separately, the conditions were 
sufficiently stable that we could use the same flux calibra- 
tion scale for the entire dataset. The zero point magni- 
tudes were K s = 22.27 ±0.02 for the February/March run 
and K s — 22.34 ± 0.03 for the October run, with uncer- 
tainties reflecting the scatter from multiple standard stars 
observed at a variety of times, airmasses, and positions on 
the detector during the course of each night. These and 
all other magnitudes in this paper are Vega-relative. 

After measuring the K s magnitudes of the central stars 
(see Table ^) , we wanted- to the extent possible- to re- 
move all traces of them from the images. To eliminate the 
(approximately) azimuthally symmetric part of each star's 
PSF, we fit ellipses to its isophotes. Inside the FWHM of 
the star's surface brightness profile, we allowed the posi- 
tion angle and ellipticity to be free parameters; outside the 
FWHM, these were fixed to the values at the FWHM to 
permit isophote-fitting to the PSF wings out to a radius 
of ~ 15". Deviating pixels (i.e., those that could in prin- 
ciple be due to faint underlying sources) were rejected in 
these fits. The azimuthally asymmetric diffraction spikes 
in the star's PSF had in general already been excised at 
the earlier stage of image combination, where their dif- 
ferent orientations for the different sets of ten frames al- 
lowed them to be eliminated via high pixel rejection. For 
the small minority of fields where this was not possible or 
that centered on very bright stars, no satisfactory method 
could be found to remove the diffraction spikes that did 
not also compromise the quality of the data. 



100 



E 
o 
u 




18 19 
K s magnitude 

Fig. 2. Completeness as a function of magnitude. The 
solid line is the mean over all 42 fields; the error bars indi- 
cate the dispersions between fields, reflecting field-to-ficld 
variations in sensitivity. 



robust against aperture effects (Vaisanen et al. 2000 
cf. Martini l2lM|l : 

its semi-major axis, defined as the second-order mo- 
ment (i.e., the intensity- weighted RMS) parallel to its 
major axis; 

its semi-minor axis, defined as the second-order mo- 
ment perpendicular to its major axis; 
the position angle of its major axis, in degrees east of 
north; and 

its SExtractor stellarity index, ranging from (diffuse) 
to 1 (pointlike). 



3. Analysis 

3.1. Object identification 

We used SExtractor (Bertin & Arnouts 1996) to generate 
object catalogues from our final, PSF-subtracted images. 
Table|211ists the sources, for each of which we provide the 
following information: 

1. a unique identifier that combines a field name from 
Table □ with the source's right ascension and decli- 
nation offsets from the central star in tenths of arc- 
seconds (SBSF 01+047-093, for example, is the ob- 
ject 4.7" east and 9.3" south of the bright star in field 
SBSF 01); 

2. its right ascension, estimated from the right ascension 
offset and columns (3) through (5) of Tabled 

3. its declination, estimated from the declination offset 
and columns (4) and (6) of Table ^ 

4. its angular separation from the bright star; 

5. its K s magnitude, using the SExtractor BEST defini- 
tion that appears in simulations to be satisfactorily 



Within Table the bright star fields are listed in order of 
right ascension; within each field, the individual sources 
are then listed in order of increasing angular separation 
from the central star, out to a maximum of A6 — 60". The 
full catalogue includes 1980 sources; this total does not 
include the 42 bright stars on which the fields themselves 
are centered. 

We estimate that the offsets of each field's sources rel- 
ative to the central star are accurate to within the size of 
a SOFI pixel, i.e., ±0.075". Although this relative astrom- 
etry formally applies only to the 2001.1 and 2001.8 epochs 
of our observations, the nominal proper motions of the ref- 
erence stars are small enough that (with the exception of 
SBSF 05) J2000.0 coordinates can essentially be computed 
from the field centers and source identifiers alone. We note 
that the absolute coordinates listed in Table [21 may be 
inaccurate if the central stars' true proper motions differ 
from the relative values listed in the USNO-B1.0 catalogue 
and Tabled or if some of the (stellar) objects themselves 
have nonzero proper motions. 
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Table 3. Sources in bright star fields identified by SExtractor; columns are described in Sect. 13.11 A full version is 
available in electronic form at the CDS. 



Source 


T> A 
ti. A. 


Dec. 


ill J 




Major axis 


Minor axis 


Ti A 

r. A. 


btellarity 


identifier 


(J2000) 


(J2000) 


(arcsec) 


(mag) 


(arcsec) 


(arcsec) 


(degrees) 


index 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


SBSF 01+077+024 


00 11 39.38 


-19 37 08.6 


8.1 


14.04 


1.63 


1.11 


81 


0.09 


SBSF 01+047-093 


00 11 39.17 


-19 37 20.3 


10.5 


15.35 


1.63 


0.65 


50 


0.03 


SBSF 01-132-076 


00 11 37.90 


-19 37 18.6 


15.2 


20.12 


0.14 


0.11 


28 


0.42 


SBSF 01+130+149 


00 11 39.76 


-19 36 56.1 


19.8 


18.60 


0.35 


0.29 


71 


0.96 


SBSF 01-217-099 


00 11 37.30 


-19 37 20.9 


23.9 


17.83 


0.49 


0.42 


131 


0.03 


SBSF 42-057-480 


23 29 55.36 


-18 36 42.0 


48.4 


19.34 


0.30 


0.22 


60 


0.55 


SBSF 42-487+056 


23 29 52.34 


-18 35 48.4 


49.0 


18.26 


0.47 


0.37 


147 


0.19 


SBSF 42+516-130 


23 29 59.40 


-18 36 07.0 


53.2 


14.12 


0.61 


0.56 


91 


0.98 


SBSF 42+335-466 


23 29 58.12 


-18 36 40.6 


57.5 


19.71 


0.23 


0.21 


170 


0.51 


SBSF 42-016+590 


23 29 55.65 


-18 34 55.0 


59.0 


20.82 


0.12 


0.10 


47 


0.38 



Table 4. Galaxy number counts as a function of K s magnitude. Column (1) indicates the magnitude bin. Column 
(2) list the corresponding completeness f c (when reliable). Columns (3), (4), and (5) list the raw numbers of galaxies 
with source separations AO < 5", 5" < AO < 30", and 30" < AO < 60", respectively. Columns (6), (7), and (8) 
list the completeness-corrected number counts of galaxies per magnitude per square degree for source separations 
5" < AO < 30", 30" < AO < 60", and 5" < A9 < 60", respectively. We do not list number counts for the two faintest 
magnitude bins due to the substantial uncertainty in the true f c at these levels. Uncertainties are calculated from 
Poissonian statistics following Gehrels ([1986 ) . 





fc 


A raw (0, 5) 


A raw (5, 30) 


A raw (30, 60) 


N (5, 30) 


N (30, 60) 


N (5, 60) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


12.5 - 13.5 


1.000 





2 


1 


zzy: -145 

449ZS 


o 6 +84 
JO -30 


821™ 


13.5 - 14.0 


1.000 





2 


1 


70+167 
'°-61 


165152° 


14.0 - 14.5 


1.000 





4 


1 


8981SS 


70+167 
' J -61 


07c + 185 

z,i o_ 119 


14.5 - 15.0 


1.000 





4 


5 


8981S8 


364I 2 « 


494+^6 


15.0 - 15.5 


1.000 





5 


5 




36413? 


549l 2 ^ 


15.5 - 16.0 


1.000 





7 


10 


1 571 + 847 


79 o+310 
/ Zo 227 


934lf 2 t 


16.0 - 16.5 


0.994 





9 


23 


2033±^ 


16841^ 


1769+|Ii 


16.5 - 17.0 


0.990 


1 


15 


28 


34021^4 


2058l«7 


2387l*f 3 


17.0 - 17.5 


0.986 





15 


36 


34151$° 


2656lgJ 


2842+g 


17.5 - 18.0 


0.970 





30 


68 


6943tg£ 


51001SS 


5552_ 660 


18.0 - 18.5 


0.958 


1 


54 


80 


12655t$gj 


6076l« 8 8 


76861™ 


18.5 - 19.0 


0.946 





62 


142 


1 471 4+ 2117 

Lit lt_ ls63 


1092l!f?{j 


11849!!!°, 


19.0 - 19.5 


0.888 





69 


188 


17445™ 


15403tE!g 


15903+ggf 


19.5 - 20.0 


0.717 





78 


235 


244231 3 ™ 


23846lJ«g 


2398711S* 


20.0 - 20.5 


0.439 





56 


173 


28638+3^4 


28672+^so 


28664l 2 °g 


20.5 - 21.0 


0.127 


1 


32 


114 


56568™° 


653091^°? 


631701^ 


21.0 - 21.5 







3 


26 








21.5 - 22.0 







1 


2 









3.2. Completeness and photometric accuracy 

We have used a single set of simulations to determine si- 
multaneously the photometric accuracy and completeness 
limit as a function of magnitude for our full source cat- 
alogue. For each field, we extracted a 10" x 10" cutout 
around every object with K s < 18 and rescaled it by a 
constant factor so that the source magnitude was shifted 
into the range 16 < K s < 21. Brighter objects with true 



K s < 15 (about half of which are the central stars them- 
selves) were dimmed to magnitudes randomly distributed 
in the range 16 < K s < 18; fainter objects with true 
15 < A" s < 18 were dimmed to magnitudes randomly 
distributed in the range 16 < K s < 21, subject to the 
constraint of > 1 mag fading. We then pasted the rescaled 
cutout back into the original image at a random position 
with separation 5" < AO < 55" from the central star. 
We did not consider test separations A9 < 5" due to the 
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Fig. 5. Galaxy surface densities as a function of separation from the bright star in our non-NVSS fields. Left panel 
shows mean for 23 fields within EIS Wide patches; right panel shows mean for 14 fields outside of EIS Wide patches. 
Only for AO < 5" does the star have an effect on the detectability of other objects. Dashed lines indicate the average 
number densities for 5" < AO < 60". The slight increase in number density at A8 < 30" for the EIS fields results from 
our having preferentially imaged fields with one or more detected /-band sources at AO < 30". 




1 7 



18 19 
K s magnitude 



20 



21 



Fig. 3. Uncertainty in our estimate of K s magnitude as a 
function of true source magnitude. As described in Sect. 
13.21 we plot the variance of the distribution of measured 
magnitudes for each set of nominal magnitudes (i.e., of 
fictitious objects pasted into simulated images). 



complex coupling of source detection with PSF subtrac- 
tion at small radii (see e.g., Fig. EJ, but avoided separa- 
tions 55" < AO < 60" merely to minimize edge effects. 
The enforcement of > 1 mag fading ensured that the noise 
would remain identical to that in the original image. Using 
the same parameters as in the original execution, we ran 
SExtractor on the modified image and attempted to re- 
cover the fictitious object within 0.5" of its nominal po- 
sition. The procedure was repeated 100 times for every 
K s < 18 object in each of the 42 fields, giving a total of 
57,800 iterations. 

Figure El and Table 0] show the completeness we de- 
rive over the range 16 < K s < 21 from the dataset as a 



* 2 




_i_ 



14 16 18 20 

K s magnitude 

Fig. 4. Differential galaxy number counts as a func- 
tion of K s magnitude. Triangles and squares indicate 
completeness-corrected counts for the disjoint sets of 
sources with separations 5" < AO < 30" and 30" < AO < 
60" from the bright stars (i.e., columns (6) and (7) in 
Table |IJl, respectively. The thick solid line shows the rela- 
tion (also completeness-corrected) when these two subsets 
are combined (column (8) in Table 0J. 



whole, i.e., treating the sum of our 42 5" < AO < 60" 
annuli as the "discrete deep field" equivalent of a single 
contiguous 131 arcmin 2 area with patchy sensitivity. We 
estimate our catalogue is 90% complete to K s = 19.2 and 
50% complete to K s — 20.2. These values include the ef- 
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fects of non-detections due to the proximity of another 
(brighter) object, typically 1-2% at K s ~ 17, and apply 
as well to extended as to point sources. The error bars in 
Fig. ^indicate the field-to-field variations in completeness, 
which result from differences in integration time, seeing, 
and source surface density. 

Because SExtractor calculates a magnitude for each of 
the fictitious sources it recovers in the 57,800 simulated 
images described above, we can estimate the uncertainty 
in our photometry from the distribution of differences be- 
tween the recovered sources' nominal and measured mag- 
nitudes (Am = m mcas — m nom ). Empirically, this distribu- 
tion is symmetric about zero (i.e., < Am >= 0) down to 
K s = 20 with large non-Gaussian wings. At fainter levels, 
< Am > begins to skew negative as sources with m moas > 
m nom become increasingly difficult for SExtractor to re- 
cover. Figure |3| plots the variance < Am 2 > — < Am > 2 
(after rejection of a few catastrophic outliers) as a func- 
tion of m nom . The implied errors are significantly larger 
than the formal uncertainties returned by SExtractor, and 
provide a more realistic indication of the true accuracy of 
our photometry. 

3.3. Galaxy counts vs. K s magnitude 

In order to evaluate the reliability of our source cata- 
logue, we have derived a magnitude-counts relation from 
our data suitable for comparison with the published liter- 
ature. We begin by using the SExtractor stellarity index 
to perform star/galaxy separation within our source cata- 
logue. In addition to the 240 objects with stellarity index 
> 0.97 that are clearly stars, we exclude an additional 
151 sources with stellarity index above a more conserva- 
tive threshold of 0.90. This leaves 1589 sources in the 42 
catalogue fields, of which 3 lie at separations A9 < 5", 448 
at 5" < A6 < 30", and 1138 at 30" < AO < 60". Division 
into these three subsets isolates the ranges of source sep- 
aration in which the simulations of Sect. 13.21 do not con- 
strain the completeness (i.e., AO < 5"), and in which par- 
ticularly special care will be needed to model AO PSF 
variation as a function of radius (i.e., 30" < A9 < 60"; 
see Steinbring et al. 12002(1 . Tabled bins the data by mag- 
nitude, and for the range over which we have been able to 
make reliable completeness corrections (12.5 < K s < 21) 
lists the differential number counts from both of the outer 
annuli and the full 5" < A8 < 60" range, normal- 
ized to units of counts per magnitude per square degree. 
The lcr error bars on the normalized counts are derived 
from Poissonian statistics according to the prescriptions 
of Gehrels itMl) . 

Figure 01 plots the completeness-corrected counts for 
the two annuli (as discrete points with error bars) and 
for the full sample (as a continuous curve). It is immedi- 
ately apparent- as from comparison of columns (6) and 
(7) in Table 0J- that our 42 fields in aggregate have a 
higher surface density of background sources at separa- 
tions 5" < AO < 30" from the star than at separa- 



* this paper 

* Bershady et al. (1998) 
ffl Djorgovski et al. (1995) 
Gardner et al. (1993) 

n Gardner et al. (1996) 
Huang et al. (1997) 

* Huang et al. (2002) 
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16 



K magnitude 



20 



Fig. 6. Differential galaxy number counts as a function 
of K s magnitude, compared to results from the literature. 
The solid line indicates the best power-law fit to our data 
(stars with error bars) for 16 < K s < 21 and 5" < AO < 
60" and has a slope a = 0.33 ± 0.01. References for other 
symbols are as indicated in the figure. 

tions 30" < AO < 60" in the brighter magnitude bins. 
This result is not surprising, given that a number of our 
fields were selected in part precisely because they did re- 
veal background sources within 30" of the central stars 
in /-band images. In particular, when we compare plots 
of mean galaxy surface density versus angular separation 
for 23 fields in EIS Wide patches (most subject to this 
bias) and for 14 fields not in EIS Wide patches or near 
NVSS sources, as in the two panels of Fig. we see a 
slight step in surface density at ~ 30" in the former plot 
that is absent from the latter. 3 Note that this effect arises 
entirely in the bright magnitude bins; at fainter magni- 
tudes (K s > 19), Table Hand Fig. H indicate that there is 
no gradient in the source density as a function of source 
separation. At the still fainter magnitudes that can only 
be reached with AO-assisted imaging of these fields, there 
should again be no dependence of source surface density 
on separation from the bright star. 

Figure shows a comparison of the completeness- 
corrected galaxy number counts calculated from the full 
5" < AO < 60" sample with comparable results from 
the literature. This figure includes a line showing the 
best power-law fit to the number counts for the range 
16 < K s < 21; we derive a slope a = d\ogN/dm = 
0.33 ± 0.01. Leaving out the K s ~ 20.75 point, which has 



3 The slight difference in mean surface densities between EIS 
(11.3 arcmin -2 ) and non-EIS (13.9 arcmin -2 ) fields may be due 
to residual contamination by faint stars in the latter, which 
tend to lie at somewhat lower Galactic latitudes. 
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a large and uncertain correction for incompleteness, gives 
a = 0.31 ± 0.01; recomputing the counts from the 37 non- 
NVSS fields only (to circumvent the plausible clustering 
of the NVSS fields' faint sources) gives a = 0.34 ± 0.01. 
In all cases, the agreement with previously published val- 
ues in the range 0.32 — 0.37 is excellent (e.g., Djorgovski 
et al. McLeod et al. Bershady et al. 

McCracken et al. I2"UUU1 Huang et al. I2001[> . Our number 
counts are also in good agreement with prior results in 
terms of normalization for K s > 16 (although brighter 
magnitude bins are affected by the selection bias discussed 
above). The overall agreement with the literature gives us 
confidence that our approach to source identification is 
sensible, and (for faint magnitudes) that our fields define 
a representative slice of the extragalactic 2.2 fim sky. 

4. Conclusions 

The agreement of the magnitude-counts relation derived 
from analysis of our "discrete deep field" with previ- 
ous derivations by other authors is not a great surprise. 
Indeed, we would expect a priori that our results should 
be no more biased by cosmic variance than those of 
groups working with contiguous fields of equivalent to- 
tal area, although the latter will be more useful for di- 
rect measurements of clustering strengths. The consis- 
tency demonstrates that our somewhat ad hoc approach 
to field selection- with some dependence on the availabil- 
ity of EIS data and/or the proximity of a bright radio 
or optical source- has not yielded a biased picture of the 
deep 2.2 /j,m sky. As a result, it is possible to proceed to 
AO observations of this set of fields with full confidence 
that any conclusions based on the aggregate properties of 
their faint near-IR sources will be statistically sound. We 
hope that the southern hemisphere user community will 
find the catalogue of 391 faint stars and 1589 galaxies we 
have provided here to be useful in fully exploiting the use 
of 8m-class telescopes at or near their diffraction limits. 
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